Energy spectra of gamma-rays, electrons and neutrinos produced at proton-proton 

interactions in the very high energy regime 
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We present new parametrisations of energy spectra of secondary particles, 7r-mesons, gamma-rays, 
electrons and neutrinos, produced in inelastic proton-proton collisions. The simple analytical ap- 
proximations based on simulations of proton-proton interactions using the public available SIBYLL 
code, provide very good accuracy for energy distributions of secondary products in the energy range 
above 100 GeV. Generally, the recommended analytical formulae deviate from the simulated dis- 
tributions within a few percent over a large range of x = Ei/E p - the fraction of energy of the 
incident proton transferred to the secondaries. Finally, we describe an approximate procedure of 
continuation of calculations towards low energies, down to the threshold of 7r-meson production. 

PACS numbers: 13.75.Cs, 13.20. Cz, 13.60.Hb, 14.60.-z 



I. INTRODUCTION 

Any reliable interpretation of an astronomical obser- 
vation requires not only high quality experimental in- 
formation concerning the spectral, temporal and spatial 
properties of radiation, but also unambiguous identifica- 
tion of the relevant radiation processes. In this regard, 
the good knowledge of characteristics of radiation mech- 
anisms is a principal issue in astrophysics, in particular 
in gamma-ray astronomy where we often face a problem 
when the same observation can be equally well explained 
by two or more radiation processes. 

Fortunately, all basic radiation processes relevant to 
high energy gamma-ray astronomy can be comprehen- 
sively studied using the methods and tools of experimen- 
tal and theoretical physics. This concerns, in particular, 
one of the most principal gamma-ray production mech- 
anisms in high-energy astrophysics - inelastic proton- 
proton interactions with subsequent decay of the sec- 
ondary 7T° and 77-mesons into gamma-rays. The decays 
of accompanying 7r ± -mesons and some other (less impor- 
tant) short-lived secondaries result in production of high 
energy neutrinos. This establishes a deep link between 
the high energy gamma-ray and neutrino astronomies. 
Finally, the secondary electrons (positrons) from p-p in- 
teractions may compete with directly accelerated elec- 
trons and thus significantly contribute to the nonthermal 
electromagnetic radiation of gamma-ray sources from ra- 
dio to hard X-rays. 

The principal role of this process in high energy astro- 
physics was recognized long ago by the pioneers of the 
field (see e.g. UlEIIl), 

in particular in the context of 
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their applications in gamma-ray (e.g. jj|) and neutrino 
(e.g. yj) astronomies. The first reliable experimental 
results obtained with the SAS-2 and COS-B gamma-ray 
satellite missions, initiated new, more detailed studies of 
7r°-production in p — p interactions |a l7j, in particular 
for the interpretation of the diffuse galactic gamma-ray 
background emission. The prospects of neutrino astron- 
omy initiated similar calculations for high energy neutri- 
nos and their links to gamma-ray astronomy |&lll llfT , [llj . 

The observations of diffuse gamma-ray emission from 
different parts of the galactic disk by EGRET revealed 
a noticeable excess of flux at energies of several GeV. 
Although this excess can be naturally explained by as- 
suming somewhat harder proton spectrum (compared to 
the locally measured flux of cosmic rays) fl^. Ilfll H^ |. 
this result recently initiated a new study of the process 
[lfll ] based on the approach of separation of the diffrac- 
tive and non-diffractive channels of interactions as well 
as incorporating violation of the Feynman scaling. 

It is remarkable that while precise calculations of 
gamma-ray spectra require quite heavy integrations over 
differential cross-sections measured at laboratory exper- 
iments, the emissivity of gamma-rays for an arbitrary 
broad and smooth, e.g. power-law, energy distribution 
of protons can be obtained, with a quite reasonable accu- 
racy, within a simple formalism (see e.g. Ref.^3|) based 
on the assumption of a constant fraction k of energy of 
the incident proton released in the secondary gamma- 
rays (see 16] for comparison of different approaches). 
This approach relies on the energy-dependent total in- 
elastic cross-section of pp interactions, '(£?), and as- 
sumes a fixed value of the parameter k w 0.17, which 
provides the best agreement with the accurate numerical 
calculations over a wide energy range of gamma-rays. 

On the other hand, in the case of sharp spectral fea- 
tures like pileups or cutoffs in the proton energy distribu- 
tion, one has to perform accurate numerical calculations 
based on simulations of inclusive cross-sections of produc- 
tion of secondary particles. Presently three well devel- 
oped codes of simulations of p-p interactions are public 
available - Pythia [HI, SIBYLL pj|, QGSJET [Of. The 
last two as well as some other models are combined in 
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the more general CORSIKA code |2J| designed for sim- 
ulations of interactions of cosmic rays with the Earth 
atmosphere. These codes are based on phenomenological 
models of p-p interactions incorporated with comprehen- 
sive experimental data obtained at particle accelerators. 

While these codes can be directly used for calculations 
of gamma-ray spectra for any distribution of primary 
protons, it is quite useful to have simple analytical pa- 
rameterizations which not only significantly reduce the 
calculation time, but also allow better understanding of 
characteristics of secondary electrons, especially for the 
distributions of parent protons with distinct spectral fea- 
tures. This concerns, for example, such an important 
question as the extraction of the shape of the proton en- 
ergy distribution in the cutoff region based on the anal- 
ysis of the observed gamma-ray spectrum. Indeed, while 
the delta-functional approach implies, by definition, sim- 
ilar spectral shapes of gamma-rays and protons (shifted 
in the energy scale by a factor of k) , in reality the energy 
spectrum of highest energy gamma-rays appears, as we 
show below, smoother than the distribution of protons 
in the corresponding (cutoff) region. Thus, the use of 
simple approximations not always can be justified, and, 
in fact, they may cause misleading astrophysical conclu- 
sions about the energy spectra of accelerated protons. 

Two different parameterizations of inclusive cross- 
sections for pion-production in proton-proton interac- 
tions has been published by Stephens and Badwar [g] and 
Blattnig et al |2j|. However, recently it was recognized 
[lfl that both parameterizations do not describe correctly 
the gamma-ray spectra in the high energy regime. The 
parameterization of Stephens and Badwar underpredicts 
the yield of high energy pions. The parameterization by 
Blattnig et al. is valid for energies below 50 GeV; above 
this energy it significantly overpredicts the pions produc- 
tion. 

In this paper we present new parameterizations for 
high energy spectra of gamma-rays, electrons and neu- 
trinos based on the simulations of proton-proton inter- 
actions using the SIBYLL code 0], and partly (only for 
distributions of 7r-mesons) the QGSJET code. We 
provide simple analytical approximations for the energy 
spectra of secondaries with accuracies generally better 
than several percent. 



convenient to present F„(x, E p ) in the form 



F n (x, E p ) = fa*{x, E p ) 



(2) 



By definition 

dN^ ee F w {x, Ep) dx = F w (x, Ep) dE^/Ep (3) 

is the number of neutral 7r-mesons per one p — p inter- 
action in the energy interval (E„ , E v + dE^). The pre- 
sentation in the form of Eq.(0) is convenient because it 
allows to estimate easily the multiplicity of 7r-mesons. 
The function <j>(x, E p ) obviously can be presented in dif- 
ferent forms. For the simulated distributions obtained 
with the QGSJET model (see Fig. the function 
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(4) 



provides a rather good fit with the parameters , a and 



r as weak functions of E p . Note that 



J F n (x, E p ) dx — Bn . 



(5) 



Also, at the threshold, x — m n /Ep, the function 
Fn should approach zero. In order to take this ef- 
fect into account, we introduce an additional term (1 — 
m 7r /(x Ep)) 1 / 2 and require that F v = at x < m v /E p . 
Then we obtain 



F w (x, E p ) = AaBirX 



a-l 



l-x c 



(1 + rx a Y 



3r 



1 - x c 



1 + rx a 



1 - 



x E r 



1/2 



(6) 



The spectra of charged pions are described by the same 
equation. 

The parameters B„, a, and k were obtained via the 
best least squares fits to the simulated events (histogram 
in Fig. nj . For the energy interval of incident protons 
0.1 — 10 3 TeV this procedure gives 



B^ = 5.58 + 0.78 L + 0.10 L 2 , 



(7) 



II. 



INCLUSIVE SPECTRA OF PIONS 



The energy spectra of secondary products of p-p in- 
teractions F 7r (x, E p ) are expressed through the total and 
inclusive cross-sections, 



E„ f d 3 a , 9 
F„(x,E p ) = -2- / -r^d 2 Pt , 

Omel J dp J 



(1) 



where x — E^jE p is the ratio of the energy of incident 
proton Ep transferred to the secondary 7r-meson. It is 



3.1 



0.89 



Bl' 2 (1 



-0.33 B, 



) 



(8) 



where L = bx{E p /l TeV). 

For x > 10" 3 (provided that E^ > 10 GeV) this ana- 
lytical approximation deviates from simulations less than 
10 percent for the entire energy range of incident protons 
from 0.1 to 10 3 TeV. In Fig.[l]the results of the QGSJET 
simulations and the fits based on Eq.jSJ are shown for the 
proton energies of 0.1 TeV and 10 3 TeV. 
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The results of numerical simulations of the energy dis- 
tribution of secondary pions obtained with the SIBYLL 
code are well described by the function 



^sibyll — —B n 



1 - x a 



(9) 



1 + rx a (1 - x a ) 
with the best fit parameters 



98 2 6 
B^ = a + 0.25, a=-^,r = — , (10) 



where 



a = 3.67 + 0.83 L + 0.075 L 2 , L = \n(E p /l TeV) . (11) 

From Eqs.© and and introducing an additional 
term (1 — m 7r /(x Ep)) 1 / 2 , one obtains 



F w (x, E p ) = AaB^x° 



l-x° 



1 + rx a {\ - x a ) 



1 r(l-2x a ) 
! — 

l-x a 1 + rx a {\ - x a ) 



1 - 



XE r . 



1/2 



(12) 



In the context of different astrophysical applications, 
the ultimate aim of this study is the analytical descrip- 
tion of gamma-rays and lepton from decays of unstable 
secondary particles produced at proton-proton interac- 
tions. In this regard, in the case of gamma-rays one has 
to consider, in addition to 7r°-mesons, also production 
and decay of 77-mesons. The analysis of the simulations 
show that the energy spectrum of 77-mesons is well de- 
scribed by the function 



F v {x, E p ) = (0.55 + 0.028 lux) 1 



xE r , 



B^ix, Bp) 



with the condition F v (x, E p ) — at x < m v /E p . 



(13) 



In Fig. [2 we present the results of simulations of en- 
ergy distributions of n- and 77-mesons obtained with the 
SIBYLL code, together with the analytical approxima- 
tions given by Eqs. 11211 and 11311 for four energies of 
incident protons: 0.1, 10, 100, and 1000 TeV. Note that 
Eq. fl2} describes the spectra of pions at x > 10~ 3 with 
an accuracy better than 10 percent over the energy range 
of protons from 0.1 TeV to 10 5 TeV. The approximation 
for the spectrum of 77-mesons is somewhat less accurate. 
However, since the contribution of 7r°-mesons to gamma- 
rays dominates over the contribution from decays of 77- 
mesons, accuracy of Eq. Ill -it is quite acceptable. 

The QGS.IET and SIBYLL codes give quite similar, 
but not identical results. The differences are reflected in 
the corresponding analytical presentations. In FigO the 
energy spectra of pions described by Eqs.(U and l(T2|l 
are shown for two energies of protons. While these two 
approximations give similar results around x ~ 0.1, de- 
viations at very small, x < 10~ 2 , and very large, x > 0.5, 
energies of secondary pions, are quite significant. Gen- 
erally, the gamma-ray spectra from astrophysical objects 
are formed at a single p-p collision (the so-called "thin 
target" scenario), therefore the most important contribu- 
tion comes from the region x = E^/Eq ~ 0.1, where these 
two approximations agree quite well with each other. On 
the other hand, for the proton distributions with distinct 
spectral features, in particular with sharp high-energy 
cutoffs, the region of x > 0.1 plays important role in 
the formation of the corresponding high-energy spectral 
tails of secondary products. Since it appears that in this 
range the SIBYLL code describes the accelerator data 
somewhat better (Sergey Ostapchenko, private commu- 
nication), in the following sections we will use the distri- 
butions of secondary hadrons obtained with the SIBYLL 
code to parameterize the energy spectra of the final prod- 
ucts of p-p interactions - gamma-rays, electrons and neu- 
trinos. 
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Figure 2: Energy spectra of it- and 77-mesons from the numerical simulations of the SIBYLL code (histograms) and from the 
analytical presentations given by Eqs. il'2l and 113(1 for four energies of primary protons: (a) 0.1 TeV, (b) 10 TeV, (c) 100 TeV, 
(d) 1000 TeV. 



III. DECAY OF PIONS 

In this section we discuss the energy spectra of de- 
cay products of neutral and charged ultrarelativistic ir- 
mesons. 



By changing the order of integration, and integrating 
over dEL, one obtains an obvious relation 



W 7 = 2 y J 7T (E^)dE 7T =2Af 7r . 



(16) 



A. 7T U — > 77 

The energy distribution of gamma-rays from the decay 
of 7r°-mesons with spectrum J^{E 7T ) reads 

00 

Ql { El )=2juE,)^, (14) 

with total number of photons 

00 00 00 

M 1 = J Q 7 (£ 7 ) dE 1 = 2 J dE 7 J ^ J^(E V ) . (15) 







The total energy of gamma-rays also coincides with the 
total energy of pions: 

00 00 
f~y = j * Ey Qj(Ej) dE^i = J E^ J^(Et^) dE^ = . (17) 



Let's consider a simple example assuming that the 
spectrum of pions has a power-law form with a low en- 
ergy cutoff at Eq: 

rfpi_/M £ ir)°i E n >E a , , 

0, E n <E , (18) 

where A, Eg and a are arbitrary constants. For a > 2 
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Figure 3: Comparison of the energy distributions of secondary 
7r- mesons given by Eq.0J (dashed curves) and En. ll2ll (solid 
curves) calculated for two energies of parent protons E p = 
1 TeV and 100 TeV. 



the number of pions and their total energy are equal 
A „ A 



Mr = 



f 



(19) 



while the energy spectrum of gamma-rays has the follow- 
ing form 

Ey > Eq , 

(20) 



Q 1 {E 1 ) 



aE* 
2A 



aE% ' 

Above the energy Eq, the gamma-ray spectrum repeats 
the shape of the pion spectrum; below Eq the gamma-ray 
spectrum is energy- independent. 

Note that in the region E > Eq the number of gamma- 
rays of given energy E is less than the number of pions 
of same energy. Indeed, 



Q- ( (E) _ 2 
ME) a 



< 1 , at E > E . (21) 



On the other hand, since the total number of gamma-rays 
by a factor of 2 exceeds the total number of pions, one 
should expect more gamma-rays at low energies. Indeed, 
the numbers of gamma-rays above and below Eq are 



Ar y {E y > Eq) 

N 1 {E 1 < Eq) 



2A 



a(a-l)E£- 
2A 



aE c n 



Although Af^(E 1 > Eq)/Mt! = Hy < 1, the sum 

N 1 {E 1 > Eq)+N 1 {E 1 < E q ) = 2 A/; . (22) 

One may say that the "deficit" of the number of high 
energy gamma-rays is compensated by the "excess" of 
low energy gamma-rays. 



B. 7T — ► ILVp 

In the rest frame of the pion, the energy of secondary 
muons and neutrinos from the ir — > [iv^ decays are 1 



n ml + mf. 

2m- ■ ' 



_ ml - ml 



2 m n 



(23) 



The momenta of secondaries in this system are 

v\=vl=El. (24) 

In the laboratory frame of coordinates (L-frame) and for 

E^ > m n 



E„ 



1 



{EkEI+p^pD « XE„, (25) 



where 



A = 1 - ml /ml = 0.427 



(26) 



Thus, the muonic neutrinos produced at the decay 
of an ultrarelativistic pion of energy E v are distributed 
as dEyfXE^ (0 < E v < XE V ). Correspondingly, if 
the energy distribution of pions is described by function 
Jtv{Eti), for neutrinos produced at the decay ir — > fxi/^ 

oo 

Qu(E v )= J ME«)£^. (27) 

E v /X 

For example, in the case of pion distribution given by 
Eo. l|T8jl . one has 



Qu(E„) = 



X a ~ 1 A 
uE% ' 
.4 



aXE« 



E v > XEq 
E v < XEq 



(28) 



The ratio of neutrinos to pions of same energy is 

Ru = = — < 1 , at E > E Q . (29) 

Note that the factor A Q_1 significantly reduces the num- 
ber of neutrinos relative to their counterpart gamma-rays 
from decays of neutral pions. 



C. \l — > ev e Vn 

The treatment of this three-particle-decay channel is 
more complex. Since to m ^> m e , in calculations we will 



1 Hereafter we assume that the speed of light c = 1. 
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neglect the mass of electrons. Then the maximum energy 
of the muon in the L-frame is 



1 



(30) 



Thus, the spectra of all particles produced at the decay 
of muons will continue up to the energy E„ (note that 
the neutrinos from the decay of charge pions continue 
to 0.427-Ett). The formalism of calculations of energy 
distributions of secondary products from decays of muons 
is described in |flll2^. 

In 7r — > fi + Vfj, decays muons are fully polarized. In this 
case the energy and angular distribution of electrons (and 
muonic neutrinos) in the rest frame of muon is described 
by the function (see e.g. |23T |1 



9e(E'J') 



1 



3-2r=F(l-2r) cos#' 



(31) 



where r = 2E' e /m fl , 9' is the angle between the momen- 
tum of electron (neutrino) and the spin of muon. The 
signs =p of the second term in Ea, l|31JI correspond to the 
decays of /i + and , respectively. 

While the distribution of muonic neutrinos is iden- 
tical to the distribution of electrons, the distribution 
g Vll {E' v ,9') of electronic neutrinos has the following form: 



g Ue (E'9') = 



6 



1 - t) 1 =F cos ( 



(32) 



where r = 2E' V / m M In the rest frame of the muon the 
maximum energy of each particle is m M /2. The functions 
given by Eqs.lpU and 1(321) are normalized as 



(33) 



g(E', 9') dE' dQ! = 1 . 



Let's denote the angle between the muon momentum 
in the rest frame of the pion and the electron momentum 
in the rest frame of muon as 9'. Since the spin of \jT 
is parallel to the momentum, 9' = 9', while for fi + we 
have 9' = ir — 9' (the spin is anti-parallel to the momen- 
tum). Therefore the energy and angular distributions of 
electrons (expressed through 9'), from ^ + and jjT decays 
have the same form: 



1 



g e {E' e ,e') = r 2 [3 + 2r(l-2r) cos0'] 



(34) 



in the Appendix. Note that the energy distribution of 
leptons from the tt + and 7r~-meson decays are identical 
because of the CP invariance of weak interactions. In 
this regard we should note that the statement of Ref.j^l 
about the "e ± asymmetry" is not correct. 

For ultrarelativistic 7r-mesons the result can be pre- 
sented in the form of functions of variable x — E/E w , 
where E is the energy of the lepton. Let's denote the 
probability of appearance of the variable x within the 
interval (x, x + dx) as dw — f(x) dx. The analytical inte- 
gration gives the following distributions for electrons and 
muonic neutrinos: 



fe{x) = f u w (x) = g Vft (x) 6(x - r)+ 
(h%(x)+h%(x))Q(r-x), 



where r — 1 — A = {m^/m^) 
3-2r 

9^ W 



9(1 -r) 2 



3- 2r 



2 - 0.573, 
(9a; 2 - 6 In a; - 4x 3 - 5) 



(36) 



(37) 



h % ^ = 9^137)1 ( 9r2 - 6 ln r 4r " 5 ) ' ( 38 ) 



^)= (1 + ^ r ~ a) x 
[9(r + x) - 4(r 2 + rx + x 2 )] 



(39) 



<d is the Heaviside function (&(x) — 1 if x > 0, and 
6(x) = for x < 0). 
For electronic neutrinos 



Ue ( x ) = 9»e ( x ) &(x-r) + 
(hi] ) (x) + h^{x))Q(r-x), (40) 

_ 2(1 - a) 
fcW "3(l-r) 2 

[6(1 -xf +r{5 + 5x-4x 2 )) + 6rlnx] , (41) 



where 



3(1 -rf 

[(l-r)(6-7r+llr 2 -4r 3 )+6rlnr] , (42) 



For the electronic neutrino we have 



9uAE'0') 



! (1-t)(1 + cos6»') 



(35) 



It is easy to obtain, after Lorentz transformations, the 
angular and energy distributions of the decay products, 
gi(E,0), in the L-frame. The integration of gi{E, 9) over 
the solid angle gives the corresponding energy distribu- 
tions of particles fi(E). These functions are derived 



3r z 

(7r 2 - 4r 3 + 7xr - Axr 2 - 2x 2 - Ax 2 r) 
The functions are normalized 



f v W (x) dx = J f Ue (a 



idx = 1 . 



(43) 



(44) 
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At x = one has / (2) (0) = 2.214 and f v (Q) = 2.367, 
while at x 1 the functions behave as 

2(3-2r)(l~x) 3 4(1 -*) 3 

/„< 2 > ( X ) = ^ ' f"e i X ) ~ 



dw/dx 



3(1 -r) 2 



1 - r 



(45) 

It is possible to obtain simple analytical presentations 
also for the so-called Z-factors: 



Zi 2 J(a)=Jx a ~ 1 f^ ) (x)dx 
o M 



4[3 - 2r - r a (3 - 2r + a - ar)} 
a 2 (l-r) 2 (a + 2)(a + 3) 

l 

Z Vi ,{a) = fx a ~ 1 fv s (x) dx — 
o 

24 [a (1 - r) - r(l - r Q )] 
a 2 (l - r) 2 (a + l)(a + 2) (a + 3) 



(46) 



(47) 



For comparison, for 7-rays and neutrinos from the tt 
\iv^ decay, one has 



and 



Z 7 = 2/a 



(48) 



(49) 



where A is defined by Eq.(J2 

In Fig. we show the functions dw/dx for different de- 
cay products in the L-frame. Note that the distributions 
of electrons and muonic neutrinos produced at decays of 
muons simply coincides. In the same figure we show also 
the distributions of photons (from the decay 7r° — > 27 
of neutral pions) and muonic neutrinos (from the decays 
7r — * [iv ^ of charged pions). The curves corresponding to 
the distribution of muonic neutrinos are marked by sym- 
bols vjp for the neutrinos from the direct decay n — * [iv^ 
and vjp for the neutrinos from the decay // 
functions are normalized, J Q dw = 1. 
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Figure 4: Energy distributions of the secondary products 
(photons, electrons, muonic and electronic neutrinos) of de- 
cays of monoenergetic ultrarelativistic neutral and charged 
pions. All distributions are normalized, J dw — 1. 



The factor 2 implies that at the decay of 7r°-meson two 
photons are produced. Similarly for the electrons and 
electronic neutrinos 

l 

Q e {E e ) = 2 j f e {x)J n (^j^, (51) 



Qu e (E Ve ) = 2|/ 1/e (a;)4^j^. (52) 
o 

Here the factor 2 takes into account the contributions 
from both ir + and ir~ pions. 

Finally, for muonic neutrinos we have 

l 

Q„„ [Eu, ) = 2 J (4a, (x) + f u p (x)) J w (^j ^ = 



D. Energy spectra of decay products for arbitrary 
energy distributions of pions 

The spectra of secondary products of decays of pions 
of arbitrary energy distributions J- k {Et T ) (the number of 
pions in the energy interval {E v , E n + dE n ) is equal 
dN n {E n ) — Jn^E-x) dE n ) are determined by integration 
over the pion energy. In particular, for gamma-rays 

l 

Q 7 (E 7 ) = 2 jj w (^x)^. (50) 
o 





The first and second terms of Eq. I|53|) describe the muonic 
neutrinos produced through the direct decay n — > fiv^ of 
the pion and the decay of the secondary muon, respec- 
tively. Note that in Ecis, l|51Jl - llo.'it we do not distinguish 
between electrons and positrons , as well as between neu- 
trinos and antineutrinos. In fact, in p — p interactions 
the number of ir + - mesons slightly exceeds the number 
of 7r~ even at energies far from the threshold. However, 



8 



N, , /N 

eptons ' 7 



- 1 1 1 r 



T 1 1 1- 




Figure 5: The ratio of the number of leptons to photons of 
same energy for power-law distributions of pions with spectral 
index a. It is assumed that 7T°, TT+ and 7r have identical 
distributions. 



0.39, respectively. Finally, note that the spectra of elec- 
trons and electronic neutrinos are very similar, thus the 
corresponding curves in FigE] practically coincide. 



IV. ENERGY SPECTRA OF PHOTONS AND 
LEPTONS PRODUCED AT P-P COLLISIONS 



For calculations of the energy spectra of gamma-rays, 
neutrinos and electrons produced at p-p interactions one 
should substitute the distribution of pions F„(x, E p ) 
given by Eq. l|l2j into Eqs.lEOIl - E3: 



' E \ dE 
ME W ) dE, -» F w (x, E p ) dx = F v [-£-, E p ) 



(57) 

The incident proton energy E p enters this equation as a 
free parameter. 



this effect is less than the accuracy of both the measure- 
ments and our analytical approximations, therefore we 
will adopt v = v and e + = er . 

It is of practical interest to compare the spectra of the 
secondary particles for power-law distributions of pions. 
Namely, let assume that the energy distributions of 7r°, 
7r~ and 7r + mesons are distributed by Ea. lll8|i with an 
arbitrary spectral index a. Then in the energy region 
E > Eo the energy spectra of photons and leptons are 
power-law with the same spectral index a. The ratio 
of leptons to photons of same energy depends only on 
a. For the muonic neutrinos, electrons, and electronic 
neutrinos these ratios are 

l 

n r 

1 ^=X a - 1 +aJ f^x^dx, (54) 
o 

i 

^=aff e (x)x a - 1 dx, (55) 
o 

i 

^ = a Jf^( x ) x a - 1 dx. (56) 
o 

The ratios given by Eas. lj54T) . ll55l) and ll56l) as func- 
tions of a are shown in Fig. Note that only in the 
case of a — 1 the ratio Meptons/-^ coincides with the ra- 
tio of the total number of leptons to gamma-rays. With 
increase of a the lepton/photon ratio decreases. For ex- 
ample, for a = 2 one has N v jN n = 0.96, N e /N y = 0.53, 
N V JN 1 = 0.52. For a = 2.5 these ratios are 0.68, 0.40, 



A. Spectra of gamma-rays 

For calculations of gamma-ray spectra one should in- 
clude the contribution also from 77-mesons. The main 
modes of decay of 77-mesons are 26]: r? — > 27 (39.4%); 
77^3vr° (32.5%); 77^ tt+tt-jt (22.6%) and 77 -> tt+tt^ 
(5%). Approximately 3.2 gamma-ray photons are pro- 
duced per decay of 77-meson; the energy transferred to all 
gamma-rays is 0.81 E v . 

The spectrum of gamma-rays produced in the direct 
7] — * 27 decay is calculated similarly to the decay 
of 7r°-mesons. For accurate calculation of gamma-ray 
spectra from the decay chains 77 — > 3ir° — » 67 and 
?] — > 7T + 7r~7r° — > 27, one has to know, strictly speaking, 
the function of energy distribution of pions g{E®) in the 
rest frame of the 77-meson, and then calculate, through 
Lorentz transformations, the distribution of pions in the 
L-frame. However, the calculations with model functions 
show that the spectrum of gamma-rays weakly depends 
on the specific form of g(E®); for different extreme forms 
of function g(E®) the results vary within 5 percent. 

The maximum energy of pions in the rest frame of 
77-meson is SJJjjj^ = (777,^ — 3mJ)/(2m 1) ), while in the 
L-frame E^ max = (E° max +p° max ) E v /m v « 0.73 E v , 
Therefore the spectrum of gamma-rays from the decay 
chain of 77-mesons breaks at E 1 — 0.73 E p . 

The spectra of gamma-rays calculated for both ir° and 
77-meson decay channels are shown in Fig. U3 The total 
spectrum of gamma-rays, based on the simulations of 
energy distributions of ir and 77-mesons by the SIBYLL 
code, can be presented in the following simple analytical 
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form 



d 

dx 



B~, 



ln(x) 
ln(x) 



1 - 



1 + feyar^T (1 - xP~t) 

1/3 \ 4 

1 — 



1 + k 7 x^(l - x^) 



ln(x) 



4/3 7 x^ 4fc 7 /? 7 x^ (1 - 2x^> ) 
1 — at^T 1 + k^x^t (1 — x^t) 



(58) 



where x = E 1 /E p . The function -F 7 (x, E p ) implies the 
number of photons in the interval (x, x + dx) per colli- 
sion. The parameters £? 7 , /3 7 , and fe 7 depend only on 
the energy of proton. The best least squares fits to the 
numerical calculations of the spectra in the energy range 
of primary protons 0.1 TeV < E p < 10 5 TeV give 



B. 

/3 7 



- 1.30 + 0.14 L + 0.011 L 
1 



1.79 + 0.11 L + 0.008 L 2 ' 
1 

0.801 + 0.049 L + 0.014 L 2 ' 



(59) 
(60) 

(61) 



where L — m(-E p /l TeV). Note that around x ~ 0.1 the 
contribution from 77-mesons is about 25 percent. In the 
most important region of x > 10 -3 Ea. 1581 describes the 
results of numerical calculations with an accuracy better 
than a few percent provided that the energy of gamma- 
rays Ey > 1 GeV. At low energies, this analytical pre- 
sentation does not provide adequate accuracy, therefore 
it cannot be used for the estimates of the total number 
of gamma-rays. 

The numerical calculations of energy spectra of 
gamma-rays from decays of ir° and rj mesons for four 
energies of incident protons, 1, 30, 300, and 3000 TeV, 
are shown in Fig. together with the analytical presen- 
tations given by Eo. 1(58)1 . The numerical calculations as 
well as the analytical presentations of the gamma-ray 
spectra shown in Fig. are based on the simulations ob- 
tained with the SIBYLL code. Recently Hillas j2^| sug- 
gested a simple parameterization for gamma-ray spec- 
trum based on simulations of p-p interactions at pro- 
ton energies of a few tens of TeV using the QGSJET 
code: F 7 (x) = 3.06 exp (-9.47 x ' 75 ) . In Fig. Owe show 
the curve based on this parameterization, together with 
gamma-ray spectra using Ea. l|58fl calculated for proton 
energies 0.1, 100 and 1000 TeV. While there is a general 
good agreement between the curves shown in Fig. es- 
pecially around x ~ 0.1, the parameterization of Hillas 
gives somewhat steeper spectra at x > 0.1 The differ- 
ence basically comes from the different interaction codes 
used for parameterizations of gamma-ray spectra (com- 
pare the curves in FigEJ. While the Hillas parameteri- 
zation (dashed line) is based on the QGSJET code, the 
calculations for gamma-ray spectra shown in Fig. by 
solid lines correspond to the SIBYLL code. 

Eqs.(H - gHJ 



B. Energy spectra of leptons 

The calculations of energy distributions of electrons 
and neutrinos from p-p have been performed numerically, 
substituting Eq.JSU into Eqs. GHJ, and JUjJ. The 
spectra of electrons from the n — > (j, decays is well 
described by the following function 

(1 + k e (\nx) 2 ) 3 , , xs - . , 

F e (x,E p )=B e ^ r no >\ (-ln(^)) 5 , (62) 



where 



B e = 



x(l + 0.3/aA) 
1 

69.5 + 2.65 L + 0.3 L 2 
1 



(0.201 + 0.062 L + 0.00042 L 2 ) 1 / 4 ' 
0.279 + 0.141 L + 0.0172 L 2 



0.3+ (2.3 + L) 



(63) 
(64) 

(65) 



Here x = E e /E^ and L = \n(E p /l TeV). 

The spectrum of muonic neutrino from the decay of 
muon, F v p) (x, E p ), is described by the same function; in 

this case x = E v /E p . 

The spectrum of muonic neutrino produced through 
the direct decays 7r — ► /1 can be described as 



F v Mx, E p ) 



B 



, Hv) 



1 - y 



0' 



y \l + k'yP'{\-yP')) 

1 4/3'/' 4k'f3'y f3 '(l-2y f3 ') 
My) " 1-yP' ~ l + /c'/'(l-/') 

where x = E Vlx /E p , y = x/0A27, 

B' = 1.75 + 0.204 L + 0.010 L 2 , 

B' = I 

A 1.67 + 0.111 L + 0.0038L 2 ' 

k' = 1.07 - 0.086 L + 0.002 L 2 . 



, (66) 

(67) 
(68) 
(69) 



The spectrum of F <i) sharply cutoffs at x = 0.427. The 
total spectrum of muonic neutrinos F v = F (i) + F < 2 ) , 
where F <i> = F e . Also, with good accuracy (less than 

5%) F Ve %F e . 

The spectra of electrons and muonic neutrinos for dif- 
ferent energies of protons (0.1, 100, and 1000 TeV) are 
shown in Figs. and El respectively. 

Although ^ -decays strongly dominate in the lepton 
production (see e.g. [jj), the contribution through de- 
cays of other unstable secondaries, first of all K-mesons, 
is not negligible. In particular, calculations based on sim- 
ulations of all relevant channels using the SIBYLL code 
show that the total yield of electrons and neutrinos ex- 
ceeds Ea. i|62)l by approximately 10 % at x ~ 0.1. At 
very small values of x the difference increases and could 
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Figure 6: Energy spectra of gamma-rays from pp interactions corresponding to four different energies of incident protons: 
(a) 1 TeV, (b) 30 TeV, (c) 300 TeV, (d) 3000 TeV. Solid lines - numerical calculations based on simulations of production of 
7r and rj mesons using the SIBYLL code, dot-dashed line - analytical presentation given by Ea. 15811 . The partial contributions 
from 7T° and r) meson decays are also shown. 



be as large as 20 % at x <C 0.1; at very large values of 
x, x > 0.5, it is reduced to 5 %. Thus, for power-law 
energy distributions of protons, Eos. ll62l) and ijfifill under- 
estimate the flux of neutrinos and electrons by a factor of 
1.1. For production of gamma-rays, the simulations per- 
formed with the SIBYLL code show that the combined 
contribution from n° and r\ meson decays exceeds 95 % 
of the total gamma-ray production rate. 



smooth (e.g. power-law type) energy distribution of pro- 
tons the so-called delta-functional approach of calculat- 
ing gamma-ray spectra gives quite accurate results (e.g. 
HiJLLOl)- However for proton spectra with distinct spec- 
tral features, like sharp pileups or cutoffs, this method 
may lead to wrong conclusions concerning, in particular, 
the predictions of the energy spectra of gamma-rays and 
neutrinos. 



To compare the contributions of the final products of 
decays, in Fig. we show the spectra of gamma-rays, 
electrons and muonic neutrinos for two energies of inci- 
dent protons, E p = 0.1 and 1000 TeV. It is interesting 
that for significantly different energies of primary protons 
the spectra of all secondary products are rather similar, 
although not identical. This explains why for broad and 
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Figure 7: Energy spectra of gamma-rays described by Ea. l58ll 
for three energies of incident protons: 0.1 TeV (curve 1), 
100 TeV (curve 2) and 1000 TeV (curve 3). The dashed curve 
corresponds to the Hillas parameterization of the spectra ob- 
tained for proton energies of several tens of TeV. 
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Figure 8: Energy spectra of electrons described by Eo. lti2fl for 
three energies of incident protons: 0.1 TeV (curve 1), 100 TeV 
(curve 2) and 1000 TeV (curve 3). 
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Figure 9: Energy spectra of all muonic neutrinos described 
by En. lti2H and Hfifijl for three energies of incident protons: 
0.1 TeV (curve 1), 100 TeV (curve 2) and 1000 TeV (curve 3). 



Let denote by 



dN p = Jp(E p )dE p 



(70) 



the number of protons in a unite volume in the energy 
interval (E p , E p + dE p ). Below for the function J p (E p ) 
we will use units cm -3 TeV -1 . Then the function 



$ 7 (£ 7 , 



d7V 7 
~dK l 



c n 



H 



n- ] , l ,, ] iE i , ) J, A E_„) F- ( -j^-, Ep 



dEp 
E„ 



(71) 



describes the gamma-ray production rate in the energy 
interval (E 7 , E 1 + dE^), where n H is the density of the 
ambient hydrogen gas, ai ne \(E p ) is the cross-section of in- 
elastic p-p interactions, and c is the speed of light. The 
function F 1 {x, E p ) is defined by En. l(58ll . For the vari- 
able x = E-y/Ep Eq.ljnj can be written in the following 
form 

$ 7 (£ 7 ) = 



V. PRODUCTION RATES AND SPECTRA OF 
PHOTONS AND LEPTONS FOR WIDE ENERGY 
DISTRIBUTIONS OF PROTONS 



The simple analytical approximations presented in the 
previous section for the energy spectra of secondary par- 
ticles produced by a proton of fixed energy make the 
calculations of production rates and spectra of gamma- 
rays and leptons for an arbitrary energy distribution of 
protons quite simple. Below we assume that the gas den- 
sity as well as the magnetic of the ambient medium are 
sufficiently low, so all secondary products decay before 
interacting with the gas and the magnetic field. 



cn H / a ine i(E 7 /x) J p {E 1 /x) F 1 (x, E 1 /x) 



dx 



(72) 



Analogous equations describe the production of neu- 
trinos v + v and electrons e - + e + ; for muonic neutri- 
nos both contributions from Ens. l(fi2|l and ijfifij) should be 
included in calculations. The inelastic part of the to- 
tal cross-section of p-p interactions o lTLQ \(Ep) can be pre- 
sented in the following form 



CTi„ci(£ P ) = 34.3+ 1.88 L + 0.25 L 2 , mb, 



(73) 



where L = ln(_E p /l TeV). This approximation is ob- 
tained with the fit of the numerical data included in the 
SIBYLL code. 
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Figure 10: Energy spectra of all decay products produced at p-p interactions for two energies of incident protons: (a) 0.1 TeV 
and (b) 1000 TeV. The spectrum of electronic neutrinos is not shown because it practically coincides with the spectrum of 
electrons. 



Thus, the calculation of production rates of gamma- 
rays, electrons and neutrinos for an arbitrary energy dis- 
tribution of protons is reduced to one-dimensional in- 
tegrals like the one given by En ^7211 . The analytical 
presentations described in the previous section can be 
used, however, only at high energies: E p > 0.1 TeV, and 
Xi = Ei/Ep > 10~ 3 . Continuation of calculations to 
lower energies requires a special treatment which is be- 
yond the scope of this paper. This energy region has been 
comprehensively studied by Dermer pf and recently by 
Kamae et al. [15]. On the other hand one may sug- 
gest a simple approach which would allow to continue 
the calculations, with a reasonable accuracy, down to 
the threshold energies of production of particles at pp 
interactions. In particular, for distributions of protons 
presented in the form 



M e p) = ex P 



(74) 



the spectra of gamma-rays and other secondaries can be 
continued to low energies u sing the ^-functional approx- 
imation as proposed in Ref.[13j. We suggest a modified 
version of this approach. Namely, we adopt for the pro- 
duction rate of 7r-mesons 



= n 6 (En - - E Mn 
n 



(75) 



where £"kin — E p — m p is the kinetic energy of protons. 
The physical meaning of parameters n and k is clear from 
the following relations: 



F^ dE-K = n , and 



kE 



kin • 



(76) 



n is the number of produced pions for the given distri- 
bution function F v , and k is the fraction of kinetic en- 
ergy of the proton transferred to gamma-rays or leptons. 



For example for power-law type distribution functions 
Fn{E) oc E~ a with a > 2, the parameter n ~ 1. 

Assuming that the parameters n and k depend weakly 
on the proton energy, from Eci J75|l one finds the produc- 
tion rate of 7r°-mesons 

, s - cn H ( E n \ ( E n 
qn{E n ) = n er incl I m p + — I J p I m p + — 

""(77) 

where K v — n/n. For the procedure described below n 
and K-K are free parameters. 

The emissivity of gamma-rays is related to q^^E^) 
though the equation 



dNi =2 f <lAEir) 
dE 1 J y/E*-mi 

£ m i n 



: dE^ , 



(78) 



where E min = E 1 + rrv^/AE^. 

The feasibility of the 5-functional approximation in the 
energy range E < 100 GeV is explained by the following 
reasons. 

1. In the energy range 1 <C E < 100 GeV the cross- 
section given by Ea J73Jl is almost constant, and the spec- 
trum of protons given by Ea. l|74fl has a power-law form. 
Therefore the spectra of gamma-rays and leptons are also 
power-law with the same index a. In this case the 5- 
functional approximation leads to power-law spectra for 
any choice of parameters n and K„. Therefore for the 
given K„ and defining the value of h from the condition 
of continuity of the spectrum at the point E = 100 GeV, 
one can obtain correct dependence and absolute value of 
the gamma-ray spectrum at 1 <C E < 100 GeV. 

2. For the value of K n = 0.17, the ^-functional approx- 
imation for power-law proton spectra agrees quite well, as 
is demonstrated in Ref . [T^| . with numerical Monte Carlo 
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calculations ^3], even at energies as low as E ~ 1 GeV 
(see also discussion in At lower energies one has to 

use, instead of Eq. I|73|) . a more accurate approximation 
for the inelastic cross-section: 



o-ineii(-Ep) = (34.3 + 1.88 L + 0.25 1? 



E 



■ mb , 



(79) 



where E th = rn p + 2m n + m 2 r /2m p = 1.22 • 10~ J TeV is 
the threshold energy of production of 7r°-mesons. Eq J79JI 
correctly describes the cross-section also at energies close 
to the threshold, and at E p > 3 E t h almost coincides with 
Ea. l|73ll . The comparison with experimental data j2^| 
shows that Eci. lf79l) can be used in wider energy range of 
protons, as it is demonstrated in Fig.lTTl 
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Figure 11: Inelastic cross-secton of p-p interactions ap- 
proximated by En. 17911 . The experimental data are from 
http:wwwppds.ihep.su:8001/c5-5A.HTML, the open points 
correspond to the cross-sections which are used in the SIBYLL 
code. 

In Fig. El we show the spectra of gamma-rays and 
leptons calculated for the proton distribution given by 
En.l|74[L The constant A is determined from the condi- 
tion 



ITcV 



E p Jp(Ep) dE p — 1 erg cm 



(80) 



In the energy range E > 0.1 TeV calculations are 
performed using Eq. lf72)l with functions Fj(x, E p ) pre- 
sented in Section ITvl at lower energies the (5-functional 
approximation is used with K„ = 0.17. As discussed 
above, n is treated as a free parameter determined from 
the condition to match the spectrum based on accurate 
calculations at E — 0.1 TeV. The parameter n depends 
on the spectrum of primary protons. For example, for 



gamma-rays, n — 1.10; 0.86, and 0.91 for power-law pro- 
ton spectra with a — 2 ; 2.5, and 3, respectively. For 
electrons these numbers are somewhat different, namely 
n = 0.77; 0.62, and 0.67 for the same spectral indices of 
protons. 

In In Fig. El we show also the gamma-ray spectrum 
obtained with the (5-functional approximation in the en- 
tire energy range (dashed line). At low energies where 
the gamma-ray spectrum behaves like power-law, the 
(5-functional approximation agrees very well with accu- 
rate calculations. However in the high energy regime 
where the proton spectrum deviates from power-law, the 
5-functional approximation fails to describe correctly the 
gamma-ray spectrum. 

It should be noted that the noticeable spectral feature 
(hardening) around E ~ 0.1 TeV in Fig. 112b is not a com- 
putational effect connected with the transition from the 
5-functional approximation at low energies to accurate 
calculations at higher energies. This behavior is caused, 
in fact, by the increase of the cross-section <7i nc i which 
becomes significant at energies of protons above 1 TeV 
(see Fig. EJ- This is demonstrated in Fig. El Indeed, 
it can be seen that under a formal assumption of energy- 
independent cross-section, erj nc i = const, the effect of the 
spectral hardening disappears. 

One can seen from Figs. IT2h .h that the exponential 
cutoff in the spectrum of protons at Eq = 10 3 TeV has 
an impact on the spectrum of gamma-rays already at 
10 TeV, i.e. as early as 0.01i?o- The reason for such an 
interesting behavior can be understood from the following 
qualitative estimates. The integrand of Eci. l|72H at £" 7 <C 
Eq has sharp maximum in the region of small x. To find 
the location of the maximum, let's assume <7i nc i = const, 
and replace the distribution F 7 from Eq. l(58l) by a model 
function Ey cx (1 — x) 4 /x which correctly describes the 
features of Ey at points x — and x — 1 . Then for a = 2 
and /? = 1 the integrand is proportional to the following 
function 



g{x) = {l-xfe-P' x . 



(81) 



where p — E 7 /Eo- Presenting the function g(x) in the 
form g(x) = exp[/i(x)], one can estimate the integral of 
this function using the saddle-point integration method, 



g(x) dx ~ exp[/i(x*)] 



(82) 



where x* = | (\/ P 2 + 16p— p) is the location of the max- 
imum of the function h(x). At p < 1 one has h(x«) ~ 
—Ap 1 / 2 — p/2. Therefore the spectrum of gamma-rays 
contains an exponential term exp [—(16 Ey/Eo) 1 / 2 ] . The 
impact of this term is noticeable already at energies 
Ey > 10~ 2 Eq. This formula also shows that the spec- 
trum of gamma-rays decreases slower compared to the 
spectrum of parent protons. In the ^-functional approxi- 
mation the gamma-ray spectrum repeats the shape of the 
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proton spectrum. The deviation of the (5-functional ap- 
proximation from the correct description of the gamma- 
ray spectrum is clearly seen in Figs lT2h .b. 

Finally we note that for a = 2 and at energies E <C Eq, 
the amplitide of the gamma-ray spectrum exceeds the 
level of the flux of muonic neutrinos. The difference is 
due to the additional contribution of 77-mesons in pro- 
duction of gamma-rays. If one takes into account only 
gamma-rays from decays of 7r°-mesons, the spectra of 
gamma-rays and muonic neutrinos are almost identical 
at energies well below the cutoff region, where both have 
power-law behavior. This conclusion is in good agree- 
ment with early studies 0, El O, |23, and does not 
contradict to the fact that in each pp interaction the 
number of muonic neutrinos is a factor of 2 larger, on 
average, than the number of gamma-rays. The disbal- 
ance is compensated, in fact, at low energies. In Fig. ITU 
we show spectra of photons and leptons at very low en- 
ergies (see Section iHlll . It is seen that the number of 
muonic neutrinos indeed significantly exceeds the num- 
ber of gamma-rays. 



VI. SUMMARY 

We present new parameterizations of energy spectra of 
gamma-rays, electrons and neutrinos produced in proton- 
proton collisions. The results are based on the spec- 
tral fits of secondaries produced in inelastic p-p interac- 
tions simulated by the SIBYLL code. The recommended 
simple analytical presentations provide accuracy better 
than 10 percent over the energy range of parent protons 
0.1 - 10 5 TeV and x = Ei / E p > 10~ 3 . 
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Appendix A: ENERGY DISTRIBUTIONS OF 
NEUTRINOS AND ELECTRONS FROM DECAYS 
OF CHARGED PIONS 

For derivation of the energy distributions of electrons 
and neutrinos at decays of charged pions, below we in- 
troduce three systems of reference - K', Kq, and K cor- 



responding to the the muon rest system, the pion rest 
system and the laboratory system, respectively. For cer- 
tainty, we will consider the 7r~-decay Since m M 3> m e we 
will neglect the mass of the electron. 

Let's denote by g e (E' e , cos 6') dE' e dCl' the energy and 
and angular distribution of electrons in the system K' , 
where 6' - is the angle between the muon spin and the 
electron momentum. Since at the ir~ decay the muon's 
spin is parallel to its momentum, for derivation of the dis- 
tribution in the K system one needs to perform Lorentz 
transformation towards the muon's spin. This results in 



g^ (E e0l cos 9 )dE e0 dfl 



1 



7^(1 - Pficos6 ) 



g e \E e ^i^(l - cos6» ), 1 C ° S ° — ) dE e0 dn , 
\ 1 — Pp. cos u J 

(Al) 

where E e o is the energy of the electron, 9o is the angle 
between the electron and muon momenta, and 



_ 1 + r _ 1 — r 

2yr 1 + r 



(A2) 



are the Lorentz factor and speed of the muon - all in the 
Kq system; r = mi/m^. 

In the Kq system the muon can be emitted, with same 
probability, at any angle, therefore Eo. ljAlll should be 
averaged over the muon directions. The function g e is 
given by Ea. lETTjl in Sec lIIII Then, denoting cos6> = 
n e n M (n e and n M - are unite vectors towards the electron 
and muon momenta) , and integrating over the directions 
of the vector n M , one obtains 

.g e0 (E e0 ) dE e0 dtl = F e (0 d£ ^ . (A3) 

47T 

Here 

F e (o = e(r - + *i 2) (0 e(£ - r ) , (A4) 

where £ = lE^jm^ varies within < £ < 1, and 



2£ 2 
37 



2 (3 + 6r - It - 4r0 , (A5) 



The spherically symmetric function given by En . (| A3 jl 
describes the energy and angular distribution of electrons 
in the pion rest system. 

Similar calculations for electronic neutrinos result in 

g^o(Eeo) dE Ve0 dn - F„ c (£) d£ ^ . (A7) 

47T 

In this equation £ = 2£ , „ e /m 7r (0 < £ < 1), and 

^.(0 = *£ } (0 6(r - + m - r) , (A8) 
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E, TeV 

Figure 13: Comparison of energy spectra of gamma-rays 
calculated for the nominal inelastic cross-section given by 
Eq.lJJ^J and a formally assumed energy independent cross- 
section Uinci = 34 mb. The proton distribution is given by 
En. 17411 . the relevant parameters are shown in figure. 



where 

= ^(3r-£-2r0, (A9) 



(2) = 4 (-r + 3£ - &e + + 3g 3 - 2rg 3 ) 
l§; (1 - r) 2 

(A10) 

Now one can calculate, for the distributions of elec- 
trons and neutrinos in the K system given by Ens JA3j) 
and (S3) the corresponding energy distributions in the 
Laboratory system K. Generally, the exact analytical 
presentations of these distributions appear quite com- 
plex, therefore in Sec lIIII we present Ecis. l^-ifijl - 14-tll de- 
rived for the most interesting case of ultrarelativistic pi- 
ons. 

Finally we should note that the energy distributions 
of electrons published in Ref.[24] are not correct. From 
derivation of Eq.(C2) of that paper one can conclude 
that this equation is obtained under the assumption that 
the muon and pion momenta are strictly parallel, i.e. 
the important step of integration of Eq.(C2) over the 
muon directions was erroneously skipped. Consequently, 
the integration of Eq.(C2) over the energy distribution 
of muons leads to inaccurate Equations (C6) - (C12) of 
Ref.Q. 
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